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TRAJECTORY CONTROL IN RENDEZVOUS PROBLEMS 
USING PROPORTIONAL NAVIGATION 
By Luigi S. Cicolani 


SUMMARY 


The rendezvous problem is defined by the end conditions that the 
position and velocity of a vehicle and its target are to be matched. In 
its present form proportional navigation theory allows the Interception 
of a target by the vehicle, that is, the matching of positions. This 
report extends the theory to include the full rendezvous end conditions. 
Trajectory constraint equations are derived and the method of computing 
the required thrust program in idealized problems is outlined. The 
thrust program obtained results entirely from the trajectory and does 
not consider the dynamics of the vehicles or errors inherent in a real 
system. The properties of the thrust program, such as its variation both 
in magnitude and direction, are examined. An acceleration forcing function 
is also derived and its properties examined. The theory is applied to 
the satellite rendezvous problem as an example and some computations are 
presented. 


INTRODUCTION 


A rendezvous between a maneuverable vehicle and its target requires 
that the position and velocity of the vehicle be matched simultaneously 
with those of the target. Such a maneuver is potentially useful in 
several commonly mentioned astronautic operations, for example, satellite 
rendezvous and planetary landings. There are also many ordinary instances 
of maneuvers requiring the rendezvous end conditions, for example, airplane 
landings, vehicle parking, refueling operations, etc. 

The control of both position and velocity to effect the rendezvous 
implies complex terminal maneuvering for which the utility of a terminal 
guidance system located in either the vehicle or target is evident. A 
number of terminal guidance systems, both manual (e.g., ref. 1 and ref- 
erences therein) and automatic (refs. 2, 3> 4), have already been proposed 
in connection with the satellite rendezvous problem. 

Reference 3 views the problem as one of orbital transfer and the 
dynamic equations of motion are solved to find the single impulse transfer 
which intercepts the target satellite. At interception a second impulse 
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reduces the relative velocity to zero. This work has been amplified in 
references 5 > 6, and elsewhere, but is essentially restricted to the 
special case in which the target is in a circular orbit. The effect of 
gravity has been linearized in this work so that the solution is in 
error if long transfer times are used, and in practice a number of 
additional corrective impulses will be required. 

Reference i)- proposes a forcing function for a continuous thrust 
automatic feedback system. This type of solution for the rendezvous 
problem has also been followed by others (ref. 7 and elsewhere). A number 
of cases have been integrated to show that the rendezvous does occur. 

But it is not clear that a rendezvous is always produced or what effect 
the external forces have on the operation of such a system, nor is it 
clear what limits must be imposed on the sensitivities of the forcing 
functions. This type of solution suffers from the analytical difficulty 
that all information on the nature of the solution can be obtained only 
by numerical integration of the dynamic equations of motion over the 
whole range of cases of interest in each rendezvous problem. Once it 
becomes necessary to integrate the dynamic equations in this problem, 
the analysis will generally yield nothing further except by an extensive 
computational program. 

There is, however, an alternative approach to the rendezvous problem 
which circumvents this analytical difficulty. It should be recognized 
at this point that the general rendezvous requirements are simply end 
conditions on the kinematics and that any practicable maneuver which 
achieves these end conditions is of possible use. The problem must 
therefore be capable of a general kinematic solution and the complete 
nature of the solution should then yield to analysis. Therefore, the 
line of investigation taken in this report is to consider first the 
kinematics and to seek a trajectory control technique which, in general, 
produces a rendezvous, leaving until later the matter of forces required 
to constrain the trajectory in the prescribed manner. 

The automatic trajectory control technique of proportional navigation 
has been applied successfully to aeronautical interception problems, that 
is, problems in which only the positions of a vehicle and its target are 
to be matched. Now the proportional navigation control equation implies 
an interception but does not of itself specify the type of interception 
any further. Thus, the homing missile maneuver of aeronautics is a special 
class of proportional navigation interceptions having the additional 
constraint that the relative velocity remain nearly constant during the 
maneuver. Such interceptions may also be used in rendezvous problems, 
provided a terminal burst of retrothrust is applied to reduce the relative 
velocity to zero. But the point here is that a rendezvous is an inter- 
ception having the constraint that the relative velocity be zero at the 
time of interception, and that there may be a general class of proportional 
navigation interceptions having this constraint. Specifically, the 
objective of this work is to find the class of proportional navigation 
interception trajectories in which the velocities of the maneuverable 
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vehicle and the target are matched, at the point of interception. The 
key to finding this class is to replace the constraint of constant 
relative velocity used in homing missile work with the velocity end 
condition of rendezvous. 

Proportional navigation is, of course, not the only trajectory 
control technique available, but the extensive experience with homing 
missile systems may be of direct benefit in rendezvous guidance systems. 

The use of proportional navigation in satellite rendezvous operations 
was originally suggested by Wrigley in referer ce 2 where the basic homing 
missile technique was applied directly to rendezvous. This involved, 
therefore, the assumption of constant relative velocity and also of small 
lead angle. Both assumptions will be eliminated; the first since, as 
explained above, we seek a more general class of rendezvous maneuvers, 
and the second because it is necessary that the initial conditions be 
unrestricted if the solution is to be of practical value. 

The properties of the thrust program required for the vehicle to 
perform the maneuvers are also investigated. Tb method for computing 
the thrust program is outlined. Once the initial conditions are 
established, the entire desired rendezvous trajectory may be computed 
and this information utilized in the dynamic equations of motion to 
compute the required thrust. The equation for thrust can be put into a 
form appropriate to a forcing function for an automatic feedback system. 
The thrust program obtained in this report, however, results entirely 
from the trajectory, and does not consider the dynamics of the vehicle 
nor the measuring and performance errors inherent in a particular real 
system. Its examination yields an insight into the general properties 
of the acceleration time histories of a successful rendezvous. 

In the last section of the report, the theory is applied as an 
example to the satellite rendezvous problem. Calculations were made for 
coplanar rendezvous with satellites in circular orbits and the character- 
istics of the thrust program are examined. Propellant requirements for 
the cases investigated are given, but a computational program s fficient 
to study fuel optimization was outside the scope of this work. 


SYMBOLS 


A 1 *,A 2 *,a 1 , 

a 2 j a 3 j b i , 

< 3-2 


ad hoc symbols in derivations 


a 


location of reference frame or origin of reference frame 


b 


observed point 



4 


f 

f T 

*V 

g s 

I 

L 

M 

Mf 

H 

R 

S,K 

t 

t f 

T 


Ul,U 2 ,U 3 




Uy 

V 

Y 
W R 
Wla 

Wy 


Xl,X 2 ,X3 

a 


4q; ^*o 


rocket exhaust velocity 

thrust force per unit mass 

target inertial acceleration 

external force per unit mass on the vehicle 

acceleration due to gravity at the satellite position 

inertial origin 

lead angle 

vehicle mass 

vehicle mass at the end of the rendezvous maneuver 

position vector of observed point 

range of observed point 

guidance system parameters 

time from start of guidance 

time of rendezvous 

guidance thrust force on the vehicle 

orthogonal unit vectors defining reference frame 

orthogonal unit vectors based on motion of observed point 

direction of relative velocity 

relative velocity vector 

relative velocity 

angular velocity of the line of sight 
inertial angular velocity of reference frame 
angular velocity of relative velocity vector 
direction cosines for the position vector, u R 
thrust angle 
Euler angles 
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7 

departure of line of sight from initial 

line of sight 

®s 

satellite orbital angular velocity 


T 

dimensionless time, 1 - ~ 

tf 


i 

transformation matrix between u R ,Ug ,u^ 
reference frame 

frame and 

% 

angular velocity of the line of sight 


% 

angular velocity of the relative velocity vector 
Subscripts and Superscripts 

( )o 

initial conditions 


f) 

vectors 


( )• 

derivative with respect to t 


o 

time derivative in reference space 


> 

time derivative in inertial space 


( )* 

quantities made dimensionless with c0 s 
THEORY 



Kinematics and Trajectory Constraint Equations 


The condition for a rendezvous is that the positions and velocities 
of a rendezvous vehicle and its destination, or target, be matched simul- 
taneously. This condition is readily shown to be independent of the 
reference frame in which velocities are measured. In particular, if the 
reference frame is based at either vehicle or target, then rendezvous 
requires that relative velocity and range become zero simultaneously, or: 

V = 0 when R = 0 

Since any trajectory of the rendezvous vehicle which has these end 
conditions will produce a rendezvous, it is apparent that one approach 
to the problem is to specify the trajectory beforehand so that the end 
conditions are met. This may be done without reference to the forces 
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involved. The thrust program necessary to follow the specified trajectory 
can be computed later from the required relative acceleration and the 
external forces. 

Proportional navigation .- Consider two vehicles^ a and b, one of 
which may be maneuvered to rendezvous with the other. A reference frame 
in which relative motion is measured is assumed to be carried by vehi- 
cle a. The motion variables may be defined after introducing an auxiliary 
set of orthogonal unit vectors referred to the observed motion of b. 



Sketch (a) 

The magnitudes and directions of the position and velocity (the 
adjective "relative" will be dropped in this section since all motion 
discussed is relative to the reference frame) of the observed point are 
defined by 


R — Rur 
V = Vuy 


These two vectors establish a plane from which an orthogonal set of unit 
vectors. Ur, ug, u$, may be defined: 

Ur = line-of-sight direction 




'UR^V 
sin L 


sin L s | Ur X Uy 1 
Uq = U<J) X Ur 

The lead angle, L, is the misalinement of the velocity vector from the 
negative line-of-sight direction. By definition it lies in the range 
0 < L < rt and hence does not indicate the direction of rotation of the 
line of sight. This information is contained in the definition of the 
vector u<j>, which is in the positive direction of the line-of-sight 
angular velocity vector. 


I 
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The angular velocities of the line-of -sight and velocity vectors 
are denoted by Wr and ¥y« The vector Wr is perpendicular to the plane 
of the line-of -sight and velocity vectors and may therefore be defined 
as 


Wr = Aru$ 

The condition that R = 0 at the end of terminal guidance may be 
met by the application of the basic proportional navigation constraint 
equation: 


Wy = SWr 


( 1 ) 


where S is a constant . Such a constraint, for appropriate values of S, 
continuously reduces the lead angle and thereby alines the relative 
velocity along the line of sight. The range, R, is also reduced con- 
currently with the lead angle, and an interception eventually occurs. 
Furthermore, the motion is now restricted to remain in a single plane 
relative to the reference frame. Physically, this may be seen by noting 
that equation ( l) requires the velocity vector to rotate only about u<j, 
which, by definition, is also the axis of rotation of the line of sight. 

Thus, the position, R, its rate of change, V, and its acceleration, V, 
are all coplanar, and once terminal guidance begins, the motion can take 
place only in the plane defined by uy and u^ at the start. This is also 
readily proved by showing that u$, which (by definition) is perpendicular 
to the plane of motion, is a constant direction in the reference space 
once proportional navigation (eq. (l)) or, indeed, any constraint of the 
form Wy = ny(t)u$ is assumed. 


Kinematic relations among the motion variables, which must be satis- 
fied by any motion, are derived from the vector identities: 


V = R 

• • • 

V = R 


( 2 ) 


where the vector differentiations are as seen in the reference frame. 
By application of the Coriolis theorem, we obtain 


Vuy — Ru^ + X R 


Vu-y + Wy XV — Ru^ 2RWj^ X Ur + Wj^ X (W-p X R) 4- Wj^ X R 


(3) 


Equation (l) may be introduced into equations ( 3 ) together with the 
geometrical relation 


uy = -(cos L)ur + (sin L)u£ (4) 

The coefficients of Ur and u^, after the indicated operations are carried 
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out, then yield the following three independent kinematic equations: 

B = -V cos L (5) 

RC2r = V sin L (6) 

-V cos L - sin L = R - Rfip 2 (7) 

Equation ( 6 ) together with the derivative of equation (5) may be intro- 
duced into equation ( 7 ) to give: 

L + (S-l)fl R = 0 (8) 


The velocity V may be eliminated between equations (5) and ( 6 ) and the 
resulting expression for Hr introduced into equation (8). Subsequent 
integration yields: 

/ \ S-l 


sin L 
sin Lq 


(9) 


The subscript zero denotes the initial value of the variable. Equation 
( 9 ) shows the functional relation between range and lead angle which 
typifies proportional navigation. It follows from ( 9 ) that S must 
exceed 1 for an interception to occur. Other functional relations among 
the variables may be obtained by substitutions among equations ( 5 ) 
through ( 9 ); for example, 

x x -(S-2) 

^ = (10) 
Vq flR 0 W 


_V_ = i_ f sin L V S-l 
V 0 i Q \sin L 0 J 


Rendezvous constraint.- Of the set of kinematic relations among the 
four variables, R, V, L, and S2 r, only three are independent, and another 
or second constraint is necessary before the motion is completely 
specified. We seek a constraint that will allow the complete motion to 
satisfy the full rendezvous end conditions within the assumption of 
proportional navigation. Expressed mathematically, the conditions on 
the choice are: 


-p-VO -pr > 
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(a) Equations ( 5 ) through (10) may not be violated and an interception 
must occur: 


S > 1 
L(t f ) = 0 

for 0 < t < t f 


dL 

at - 0 


R _ V _ L _ fiR _ 

Ro Lq £1r 


at t = 0 


where tf is defined as the time at which the range becomes zero. 

(b) A rendezvous must occur: 

V(t f ) = 0 , t f > 0 

(V(t -* 0 + ) = Vf > 0, implying an impulse at interception, is not excluded) 

These are minimum restrictions and others could be added, provided they 
are consistent with the above set. In particular, no restraints due to 
thrust limitations have been imposed. These limitations generally cannot 
be expressed as simply as the kinematic conditions. However, the con- 
ditions above are always practical rendezvous maneuvers in the sense that 
the maneuver will be completed with a finite mass. 

The second equation of constraint may be chosen variously, for 
example, forms giving any variable as a function of time, or functional 
relations among two or more variables or their derivatives. However, 
functional relations among a set of variables for which a relation is 
already specified by the preceding kinematics cannot be used since a 
new relation of this form cannot be independent without being 
c ontradi ct ory . 

It should be noted that an investigation of only one of these forms 
will suffice to cover the motion produced by a second constraint chosen 
in any other form since the choice of a second equation of constraint 
will completely specify the motion. Our main interest is to satisfy the 
requirement V(tf) = 0 so it seems natural to study velocity constraints. 
However, it is much simpler to deal with the lead angle because all 
motion variables may be given directly in terms of the lead angle and 
its derivative. Thus, a condition that is to be applied to any other 
variable has a simple equivalent condition on the lead angle or its 
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derivative. The following form is to be studied: 


where 


~ = f(r) 
■HO 



(ID 


The restrictions in (a) and (b) above will delineate the allowable 
class of constraining functions of this form. Applying these restrictions 
to f ( t) through equations (8), (9 ) , and (10) yields: 


f(0) = 0 

f ' (t) > 0 for 0 < t < 1 

f(l) = 1 

f'(l) = (S-l) t f X °- sl . n .i£ 


lim 
r -*■ 0 


ft .Li)_ = o 

r f(T )P 


-\ 




( 12 ) 


(lim 


f r (r) 


f ( t) 


— =. Vf * > 0 plus an impulse at t = 0 is also allowable) 


S-l 


where the prime indicates differentiation with respect to t. 


These conditions describe a simple class of functions. The value 
of f ( t) is prescribed at the end points, t = 0 and t = 1, and its 
derivative is specified at t = 1. Further, f r (T) is never negative 
between the end points. An additional restriction must be imposed on 
grounds that step changes in range are physically impossible. Therefore, 
step discontinuities in f(-r) are proscribed. Any function in this 
restricted class which also satisfies the last condition above will be 
satisfactory. 


The simplest function that suggests itself is 

f ( T ) = T m ( 13 ) 

(See appendix A.) So that the continuity of the present theory from 
the standard proportional navigation may be clearly exposed, a new 
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parameter, K, is defined to replace the exponent, m, by the expression 
K(S-l). The tvo parameters, S and K, become the fundamental constraint 
parameters. Applying the restrictions of (12) to equation (13) yields 
as the second equation of constraint: 

L_ = ^(S-l) 
lo 


k 

4 

9 

4 


S > 1 
K> 1 


From (12) there also follows: 


t f 


_ 5s. -4 q 

V Q sin L 0 


(14) 


(15) 


The complete motion may now be derived directly from equations (8), 
and ( 10) . 

_R_ _ f sin L Y" 1 
Bq V s ln W 


(9), 

( 16 ) 


is. = t k(s-i)-i (17) 

n R 0 

^ = ( fe)' (S ' 2V(s ' 1) - 1 <*> 

The angle 7 will denote the angle between the line of sight u r (t) and 
the initial line-of-sight direction u Ro (see sketch (c)). It is readily 
determined by noting that 


y - % 

which, from equation (17), gives: 

y = JS. Tl - tK(S-D 
S-l 


(19) 


Characteristics of the rendezvous maneuvers .- Equation ( l4) , as 
required, specifies that the lead angle decrease monotonically toward 
zero during the rendezvous maneuver. From equation (l6) the range may 
therefore increase beyond its initial value before decreasing toward zero 
if the initial lead angle is greater than jt/2, but it never increases 
if the initial lead angle is less than jt/ 2. The limiting behavior of 
the range near rendezvous (t -► 0 and therefore L is sufficiently small 
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to assume sin L = L) is readily found to be 

1 . 

R_ ~ ( . Io 
E 0 V s in Lo, 


iS-1 K 

T 


which decreases at least linearly with t. (We note that t is really 
the dimensionless time to go.) 

The angular velocity of the line of sight ( line-of -sight rate) 
given by (17) exhibits varying behavior according to the value of the 
parameter, S. In particular 


%l(tf) = 0 for S > 


K+l 

K 


2j^(tf) -*■ 00 for 1 < S < 


K+l 

K 


While motion having the latter characteristic is of doubtful practical 
interest, a valid rendezvous is still achieved. It will be seen in the 
next paragraph that the total line-of -sight angular change remains finite 
for these values of S but indefinite spiraling motion is approached 
as S -*• 1. 

Equation (19) indicates how far the line of sight moves from its 
initial position Uf> Q , as observed in the reference frame. The maximum 
departure is 

L 0 

^max g_q 


and occurs at rendezvous. This angle increases as S decreases, becoming 
arbitrarily large for S -*• 1 where an indefinite spiraling motion is 
approached. For values of S > 2 the maximum departure is always less 
than the initial lead angle. 1 

1 While off hand "loose” trajectories (S near l) do not appear to be 
desirable, these may present some advantages in a special case. Suppose 
it is desired to rendezvous with a point on a surface so that the maneuver 
ends tangent to the surface] that is, a landing. If the motion of the 
approaching vehicle is in the plane perpendicular to the landing surface, 
with the initial line of sight at an angle Aq above the surface and the 
relative velocity vector at an angle Lq below the line of sight , then 
the required value of S for a tangential landing is given from 

y = A = ifi- 
' max ^0 g_q 


where values of S < 2 are required if Lq < Aq, A vertical landing may 
be treated similarly. 
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The character of the velocity time history is obscure in 
equation (l8), which contains two uncombinable terms that may in general 
have opposing or concurrent trends depending on the values of L 0 , S, 
and K. However, some information may be obtained from the relative 
acceleration, which, after some algebra, becomes: 


V 


V 2 

-cos 


L 




tan L 
L 


Examination of the terms in this equation reveals that velocity always 
decreases during guidance if S > 2. For S < 2, velocity increases if 
the lead angle exceeds a certain value, L x , and decreases below L x . 

The value of L x depends on S and K and is given by: 


tan L x 
la 


S-2 


a K+l 
S - _ 


The limiting behavior of velocity near rendezvous (from eq. (l8)) and 
its derivative are: 

S-2 

_V_ ~ ( sin Ip NS-l K-l 

V 0 V Lo J 

2S-3 

y ~ K-l Vof_ ( sin Lp\ s-1 ^-2 

K R 0 V Lo ) 


Near rendezvous the behavior of these quantities is therefore dependent 
on K, as illustrated in sketch (b): 




r 


Velocity approaches zero continuously for values of K > 1, and at K = 1 
a step decrease of velocity occurs at t = 0. It may be anticipated 
that the thrust requirements near rendezvous will exhibit the same 
behavior as the relative acceleration, V, whose behavior depends on K. 

In particular, if K > 2, V decreases continuously to zero, and 

if K = 2, V exhibits a step decrease to zero. In the range 1 < K < 2, V 
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becomes arbitrarily large at rendezvous, and for K = 1, an impulse 
occurs corresponding to the step decrease in velocity. In the cases 
with singular behavior it can be shown that the mass ratio across the 
singular region is not zero since the integral of V remains finite 
over this region. 

For a given set of initial conditions (R 0 , V q , L d ) , the maneuver 
time, tf, (eq. ( 15 )) varies directly with K but is independent of S. 
Maneuver time is therefore a minimum for K = 1, which requires impulsive 
retrothrust at rendezvous. For K> 2 the thrust requirements near 
rendezvous are continuous, but at least twice as much time is required 
for this type of maneuver as is necessary using impulsive terminal 
thrusting. 

Equation ( 15 ) also indicates that tf becomes arbitrarily large 
for L 0 -* Jt (guidance begins with the vehicle receding along the line of 
sight from the target) or for V Q -» 0. The occurrence of these unfavor- 
able initial conditions may be obviated where necessary by providing the 
guidance system with a means of choosing the most favorable point on the 
unguided approach trajectory at which to initiate the rendezvous maneuver. 
However, the choice of initial point should consider possible fuel 
optimization as well as time requirements. 

If K = 1, then equation (l4) is the same as the corresponding 
result of reference 2 and, as must be the case, the complete motion is 
that of reference 2 generalized to include large values of lead angle. 
Thus, for this value of K the present theory reduces to the standard 
proportional navigation theory up to the point of interception. At 
interception the two differ in that the rendezvous maneuver requires a 
step decrease of velocity to zero. 


The rendezvous maneuvers derived in this section are illustrated 
in figures 1, 2, and 3- Figure 1 is a dimensionless polar plot of paths 
in the plane of relative motion (r/r 0 vs. 7 ). The origin is the origin 
of the reference frame and the vertical radial line is the initial 
line -of -sight direction. The equation governing these paths is found 
from substitutions among equations (l4), (l6), and ( 19 ): 


R 


sin[L 0 - (S-l)y] 
sin L 0 


1 

S-l 


The paths therefore depend only on Lq and S. Since the maneuver time 
is independent of S, a set of values of S will give a set of paths of 
constant maneuver time for given values of L 0 , R 0 , V Q , and K. (Changes 
in K do not affect the path for a given value of S but change only 
the time required to traverse the path.) 

Polar plots of V/V 0 versus L for various values of S and the 
initial condition L 0 = jt/2 are shown in figure 2. Here, the equation 
may be found from equations ( l4) and ( l8) . 


A 

4 

9 

4 


-f=-VO -|=- > 


2E 


15 


_V_ 

V 0 


sin L 
sin Lq. 



iy+x 

K 


f 5 - 1 


Plots of lead angle versus velocity for several values of K are given 
in figure 3. 


Reference frame «- The command trajectories n aw have been largely 
determined except to indicate the necessary transformations of vectors 
given in the the Up, ug, u^ frame to their conponents in the reference 
space of the command system. A reference space is indicated in sketch (c) 
by the orthogonal unit vectors u x , u 2 , u 3 . At this point it is not 
necessary to specify this reference frame further; it could be, for 
example, an inertial or a local vertical frame. The plane of motion of 
the observed point and its initial line of sight are established in the 
reference space by means of the Euler angles Aq, p Q , T 0 , which neces- 
sarily have constant values during guidance. The angle, 7, indicates 
the departure of the line of sight from its initial direction, uj^, as 
observed in the reference space. 


To eliminate ambiguity, Ao is to 
be measured from the positive u 3 
axis to the nearest branch of the line 
of nodes in a right-hand sense about 
ui. Similarly r o is measured from 
the same branch of the line of nodes 
in a right-hand sense about u^. 

The commuter coordinate vectors, 
up (line of sight), uq, u$, may be 
given in the satellite reference space 
by the transformation: 


( 20 ) 


“ 1 


“ “ 

% 


Ul 

% 

= V 

u 2 

_ 


u 3 

_ 


Plan# of raloflva motion 
_ (r) 



Lina of nodas 


Sketch (c) 


* = 


sin P 0 sin(r 0 +7) 
sin g 0 cos(r 0 +7) 
cos p 0 


-[sin Aocos(r 0 +7) + cos Aosin(r 0 +7)cos Pol 
sin A 0 sin(r 0 +7) - cos Aocos(r 0 +7)cos P 0 
cos Aosin P 0 


cos Aocos(r 0 +7) - sin Aosin(r 0 +7)eos f ) 0 
-[cos A 0 sin(r 0 +7) + sin A 0 cos(r 0 +7)cos g Q ] 
sin Acsln P Q 
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In general, the components of a vector in the reference space may 
be given from its components in the motion frame by the transformation: 



”Ai 


a r 

A = 

a 2 

,T 

= t 

% 


_ As _ 




T 

where \|r indicates the transpose of \| /. It will be useful to introduce 
the notation 

U R = X X U X + X 2 U 2 + X3U3 

where the Xq are the first row of ty, and to note that ug = (fip) 1 u-^ 
or 

Ug = (%)" 1 (x 1 U 1 + X 2 U 2 + X3U3) 

where the are the second row of \|/. 


A 
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Maneuver Forces 


The trajectory constraint investigation thus far has not required 
the specification of forces involved and is generally applicable to any 
problem requiring the rendezvous end conditions. However, before a 
maneuver can be carried out in any given problem, the required thrust 
program must be computed. The external forces on both vehicle and target 
and the inertial rotational motion of the reference frame must be known. 


Maneuver accelerations .- The inertial acceleration required during 
the maneuver may be found from Newton’s law. 



The reference frame is at a which may 
be either target or vehicle. If a is the 
target, then applying Newton's law to the 
vehicle we obtain: 

A + f = Jr (ha + R) 


f = 


d 2 R 

dt 2 


+ (fT - fv) 


Inertial origin 

Sketch ( d) 


or 


( 21 ) 
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f acceleration of the vehicle due to thrusting or thrust force per 
unit vehicle mass (to be referred to as the thrust program) 2 


*v 


external force per unit mass on the vehicle (gravity, aerodynamic, 
etc.) 


acceleration of the target, 


a. 2 %a 

dt 2 


The inertial relative acceleration, d 2 R/dt 2 , may be given in terms 
of the motion observed in the reference frame by the application of 
Coriolis’ theorem. Equation (2l) becomes: 

f = ? + 2W Ia X V + W Ia X (W Ia X R) + % a X R + (f T - f v ) (22) 


where _ Wj a is the inertial angular velocity of the reference frame 
and ( ) indicates time derivatives as seen in the reference frame . Since 
the desired motion is one of the rendezvous maneuvers, then the appropriate 
expressions for 7, and R from the previous text are to be introduced 
in (22), together with Wj a and the external forces. The required 
thrust program is then confuted. 

The commanded relative acceleration, V, for the rendezvous maneuvers 
may be written in terms of the measured motion variables: 


vi 

R 


(S-2)cos L - 



sin L 
1 


Uy + SWp X V 


( 23 ) 


Utilizing this form in equation (22) yields an expression for f which 
is appropriate as the basic forcing function in automatic feedback 
systems since that portion of the required thrust which depends on the 
relative motion may now be computed directly from the observed current 
motion. 


Characteristics of the maneuver acceleration .- The relative 
acceleration of the maneuvers is seen from equation (23) to consist of 
two terms. The first, directed along Uy, is largely responsible for 
satisfying the velocity end condition. The second term, SWr X V, is the 
usual proportional navigation control function which forces an inter- 
ception. Mathematically, f is apparently well behaved except possibly 
near rendezvous (t -*• 0) where both V and R become zero, and may 

sometimes be singular. Utilizing the previous results approximated for 
small values of t gives the limiting expression for V near rendezvous: 

2 If a is the vehicle, the result is: 


d g R 

dt 2 


+ (fT - fv) 


f = 


l8 


2S-3 

♦ ■ l (- ¥ ^ ♦ SLo^- 2 at X % ) (24) 

The first term, due to V, varies as (K-1 )t^-~ 2 and was previously 
discussed; it reduces continuously to zero for K > 2, becomes constant 
for K = 2, and is singular in the range 1<K<2. At K = 1 it is 
zero except at t = 0 where it becomes an impulse function corresponding 
to a step decrease in velocity of 

av = v 0 

Wear rendezvous, Uy is approximately along the negative line-of -sight 
direction so that this part of the limiting thrust program is directed 
along ur . 

The second term, due to 2 rV and directed nearly perpendicular 
to up along u$ X uy, varies as therefore reduces to zero 

continuously for S > 2/K and is singular for lower values of S. The 
singular behavior is possible within the restriction S > 1 only if K 
is less than 2. The singular behavior of the term is due to Hr modi- 
fied by the behavior of V which tends to weaken the singularity and 
reduce the range of values of S over which it may occur. In fact, 
for K > 2 the singular influence of 2 r on this term is removed. 

Because of the restriction K > 1, the maximum value of S for 
which the second term in (24) may be singular is 2. Furthermore, a 
number of possible peculiarities of the maneuver motion were noted in 
the previous text for values of S less than 2. A valid rendezvous 
is not precluded by any of these peculiarities but there is apparently 
little advantage to be derived from them. The remainder of this report 
will limit consideration to the main body of proportional navigation 
maneuvers having fairly uniform characteristics by assuming the restriction 

S > 2 (25) 

As a consequence, the second term in (24) will always reduce to zero 
more rapidly than t 2(K-1) near rendezvous. 

If Wia and frp - f v are assumed to be well behaved, as must 
generally be true in practical cases, then the first term in V' will 
be the dominant term in equation (22) near rendezvous. For K near 1, 
this dominance is delayed to values of t very close to zero because 
of the presence of (K-l) in this term; that is, the system approximates 
the constant relative velocity behavior. In these cases the behavior 
prior to the time at which the first term is dominant will be given by 
some other term in equation (22), possibly the term Wj a XV (if Wi a 4 °) 
or by frp - fy. 



-fr^O -fr" > 
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Propellant requirements .- The mass of the vehicle will vary during 
terminal guidance because of the expenditure of propellant in providing 
thrust. If the thrust system is assumed to consist of a single-stage 
rocket system of exhaust velocity c, then from Newton’s second law the 
mass ratio at any time, is: 



(26) 


where Mq is the initial vehicle mass and only the rendezvous maneuver 
forces are considered. The actual thrust at any time is then given by: 

T = M(t)f(t) (27) 


At this point, outside of the context of a specific problem, it is only 
possible to indicate those factors which affect propellant requirements 
and which therefore might be controlled so as to reduce the required 

rtf _ 

fuel. The mass ratio over the maneuver depends on the integral / |f |dt 

'-'o 

which, from (22) and (23), depends on the system parameters, initial 
conditions and reference frame . If all other parameters are specified 
in a given problem then the system parameters, S and K, could be chosen 
so as to optimize fuel consumption. However, the latitude of choice of 
these parameters will likely be limited by other considerations; for 
example, maneuver time, trajectory constraint, and thrust program 
characteristics. As previously mentioned, it is possible to choose the 
point on the approach trajectory at which guidance is begun such that the 
fuel requirements are a minimum. 


In equation (22) there are several terms which depend on the choice 
of reference frame, that is, Wp a . If the rendezvous maneuver is carried 
out in two different reference frames for an otherwise completely specified 
problem, then the two maneuvers as viewed from inertial space are 

r* f 

different. It then follows that |f| and / |f|dt will, in general, 

'“O 

differ for the two frames and that the fuel requirements will vary 
accordingly, some reference frames requiring less fuel than others in 
each particular case. While Wj a has been assumed arbitrary thus far, 
the number of useful frames is limited in practice, for example, inertial 
and local vertical frames. Of course, guidance may be carried out in 
any imaginary frame from measurements taken in one of the simpler frames 
by means of appropriate transformation relations, so that in theory an 
arbitrary Wj a is still available. For our purpose, such a procedure 
is unnecessarily complicated. Moreover, there is no simple means of 
finding the optimum Wi a . Nevertheless, for a reasonable problem 
(restricted range of expected initial conditions) it should be feasible 
with direct computations to determine if there is any significant advantage 
in the choice of Wj a from among a few simple frames. 
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As a final point, it should be noted that the thrust requirements 
for the proposed rendezvous maneuvers are, in general, variable both in 
magnitude and direction. This is illustrated in the example that follows 
dealing with the satellite rendezvous problem. 


THE SATELLITE RENDEZVOUS PROBLEM 
Application of the Force Equation 


In what follows the "target" is a satellite in planetary orbit and 
contains the command guidance system, including the reference frame. 

The reference frame is taken to be a local vertical system with Ui 
along the outward local vertical and u 3 in the positive direction of 
the satellite orbital angular velocity, that is, perpendicular to the 
orbital plane. The angular velocity of the reference frame, Wj a , is then 
the orbital angular velocity of the satellite. 

The directions of all vectors in the thrust equation (22), 
except frp - fy, are already given. The following vector relations 
are repeated here for convenience: 

R = Ru R 

V - Vuy 

Uy = -cos Lup + sin Lug 
f = Vuy + Sfi R Vng 

and also 

Wja = 6 s a 3 

The transformation relations of equation (20) now allow all terms of (22) 
to be expressed in the reference frame, giving the components of f 
along u x , u 2 , u 3 as: 


+ (^t ” ( 28 ) 
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a l 

11 

CQ 

£0 

% 

sin L 

+ V cos L 

a 2 

= R0 S - 

■ 2 0 S V 

cos L 

a 3 

II 

CQ 

£0 

% 

sin L 

+ V cos L 

bi 

II 

<! 

£0 

% 

os L - 

- V sin L 

b 2 

II 

sin L 



It remains to express the external forces, fpf - fv, in the reference 
frame. For a satellite orbit outside the planetary atmosphere, the 
difference in external forces is a result of the gravitational field of 
the planet. In a central force field, if the vehicles are sufficiently 
close during guidance, the gravity difference term may be linearized as 
shown in appendix Bj whence 



-2 0 0 


Xl 

- R 




f T ~ f V - Ss R^ 

0 10 


x 2 


_° 0 1_ 


3d. 


( 29 ) 


where R s is the radial distance of the satellite from the center of 
the force field and g s is the gravitational acceleration at the 
position of the satellite. The complete time histories of R, V, L, 2r, 
and the xp are computed from equations (lA) through (20) once the 
initial conditions and the system parameters, S and K, are given. If 
the orbit and g s are also given, then the complete thrust program may 
be computed. 


Calculations were carried out for the simple case in which the 
initial motion, and hence the entire maneuver, is coplanar with the 
orbital plane of a satellite in circular orbit around the Earth. When 
equations ( 28 ) and ( 29 ) are specialized to this simple case, then: 


0 


s 


= 0 



( 30 ) 


x 3 = x 3 = 0 


For coplanar motion the same initial position may be represented by 
two sets of Euler angles which correspond to opposite directions of 
rotation of the line of sight. These are 

Case a: Aq = ir/2 , P 0 = m/2 , f 0 

Case b: A 0 = ir/2 , P Q = -J f/2 , -F 0 
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The rotation of the line of sight is, respectively, clockwise and 
counterclockwise about the satellite axis of orbital rotation, u 3 . The 
difference between the two cases depends on the orientation of the 
relative velocity vectors (sketch (e)). 


V 


Vehicle 



(a) 



Sketch (e) 

For cases (a) and (b) the position during guidance is given by: 

xi = sin[T 0 + (±)rl 
x 2 = ~cos[r o + (±)r] 

= (±)us 


(31) 


The factor (±), which henceforth will be denoted by n, is to be taken 
as (+D or (-l) for cases (a) and (b), respectively. This procedure 
allows both cases to be included at once in the following derivation. 


Equations (2j) , (29), ( 30 ), and ( 31 ) maybe introduced in ( 28 ) to 
form an expression for the thrust program. This expression is made 
dimensionless with the factor c0 s which depends only on the rocket 
system and the particular circular orbit. The calculations will then 
apply to all circular orbits and guidance rockets. In general, these 
computational benefits are not obtained if the orbit is elliptical or 
if the gravity difference term is not linearized. 


where 

A x * 


fl* 


At* oRo^s R 

1 ^ 2 Ro 

a 2 * 


sin(r D + ny) 

f 2 * 


a 2 * 

-Ax* 


cos(r 0 + nr) 

» 4 


fs* = 0 


(32) 
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A 2 * = 



c 


, V % 

2 sm Lq — — = S- 
Vo 


cos L + 



K+l\ sin 2 Lo _R_ 1 

K / Lq Rq t_ 


Quantities made dimensionless vith c0 s are denoted by the asterisk. 

The parameters Rq0 s /c, r Q , V q /c, n, and L 0 completely specify the 
initial conditions in dimensionless form. Computations may be carried 
out in any given problem by substituting the time histories of the motion 
variables in equation (32) from equations (l4) through (19). 


Equations (26) and (27) become: 


M(t) 

-Mo” 


exp { -0 s tf / [f*]dT 


j = M(t )_ ^ 
cQgMo Mo 


(33) 


The factor 0 s t-£> is the geocentric orbit angle passed through by the 
satellite during the maneuver, or, alternatively, a dimensionless 
maneuver time. It is readily given in terms of the dimensionless initial 
conditions from equation ( 15) : 


e s t f = 



^o 

sin L 0 


Computations 


The computations covered the following ranges of initial conditions: 

n = ±1 

L q = 0 to 3rtA 

^ - 0.01 to 0.1 
c 

r Q = 0 to 2 * 

Mis = 0.OO5 to 0.05 
c 

The range of values of R<j0 s /c corresponds to initial ranges from 10 
to 100 miles for near-earth orbits (0 S = 10 -3 radians/sec) and 
c = 10,000 fps . The following values of the system parameters were taken: 



S = 5 
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K = 1, 2, 3 

These values of S and K represent fairly "tight” maneuvers (/ max = L 0 /4) 
and will illustrate the main effects of K on the thrust program. Higher 
values of K are considered unlikely because of the corresponding 
increases in maneuver time. 

The computed mass ratios for a number of sets of initial conditions 
and system parameters are given in table I. This table is intended to 
give typical values. Questions concerning values of S and K for 
minimum fuel requirements (the propellant ratio is given by 
Mp/M 0 = 1 - (Mf/Mo)) demand a more detailed computational program outside 
the scope of the present work. It appears generally from the table that 
if the fuel requirements are severe, then the variation in fuel require- 
ments with S and K will be large ~13 percent over the range of S and K 
covered for the case (R^g/c = 0.05, Tq/c = O.Ol) of table l(a). In 
such cases consideration must be given to minimizing fuel requirements 
by the proper choice of S and K. 

Two cases were computed to compare fuel requirements using the 
proportional navigation maneuvers and the two -impulse method of 
reference 3- The mass ratios, for the two methods are given in 

table II. These two isolated cases indicate comparable fuel requirements 
for the two methods, but it would be hazardous to extend this conclusion 
to apply generally. 

In figures 4 through 6, a number of thrust programs (magnitude and 
direction) are plotted versus the dimensionless time, r. The thrust 
per unit mass is given in the dimensionless form f/c6 s . Noting that 
for near-earth orbits 0 S = 10 -3 radian/sec and taking c = 10,000 fps, 
then the order of magnitude of the required acceleration, f, in feet 
per sec 2 , is given by multiplying the value of f* in the graphs by 10. 
The thrust direction is denoted by a which is the angle between the 
negative line-of-sight direction and the thrust vector, as shown in 
sketch (f) 



Sketch (f) 


a = 


r o 


+ ny + arc tan 


fi* 


(34) 
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Before discussing the thrust equation and the results, it is 
convenient to rearrange equation (22), utilizing (23), (29), (30), 
and (31) as before, except that it is now helpful to make the result 
dimensionless with V o 0 * s rather than c0 s . The result is then: 


f 

V O 0 S 


Vq 
BqQ s 



cos L - 




sin l! - 

“T~J u v 


+ S S in Lo h ^ nfl3 x 5 V +2 ^ 5s x ^ 

- 3 ^ ~ sin (r Q + n 7 )u x ( 35 ) 

*0 Ro 

The first term is due to V and is directed along uy. The second 
and third terms are directed perpendicular to uy in the same direction 
if n = +1 and in opposite directions if n is - 1 . These two terms are 
due, respectively, to SWp> X V and 2 Wj a X V. The first three terms are 
independent of the line-of -sight direction (r o + n 7 ) and therefore for 
cases in which these terms are dominant there will be very little 
variation of the thrust program with different directions of the line 
of sight . 

The fourth term is the combined effect of Wj a X (Wj a X R) and the 
gravity difference and is the only term to vary with T 0 . The importance 
of this term depends on the value of the dimensionless parameter ^oMc^s> 
if this parameter is sufficiently large (large initial velocity, low 
initial range), then the first three terms are dominant and the fourth 
term is negligible. If the parameter is small (low V Q , large Ro)? 
then the fourth term contributes important position-dependent effects 
to the thrust program. Equivalently, if Vo/RcPs is large, then the 
required maneuver acceleration is large compared to the gravity 
difference, while low V 0 /Rc4g indicates that the required acceleration 
is very low and of the same order as the gravity difference. The effect 
of Vq/R^s on the relative importance of terms is complicated by 
effects of L 0 , n, and K, but the general significance of Vo/R^'g is 
as stated above. 

Table III illustrates a case of large V 0 /Rq 9 s for two different 
positions, r 0 . A comparison of the two shows that the mass ratio and 
both the magnitude and the direction of thrust are very little affected 
by the difference in r o . Figure 4 (a) illustrates a case with low Vo/RqQs 
for two position angles. Strong effects of position and very low values 
of acceleration (about 0.6 and 0.1 ft/sec 2 for the two cases at t = l) 
should be noted. An intermediate case is given in figure 4 (b). 

Figure 5 (a) illustrates the effects of K. In the thrust program 
the value of K affects primarily the importance of the term Vuy at 
low lead angles. This effect is always seen for low values of t since 
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the lead angle will have been reduced to small values ; that Is, it is 
always seen in the limiting behavior of the thrust program. The case 
plotted in figure 5( a ) is one with a low initial lead angle (10°) so that 
the influence of K is pronounced throughout _most of the maneuver. 

For K = 1, the chief function of the term Vuy is to remove the 
relative velocity at t = 0 by means of an impulse, noted by a break in 
the f* plot of figure 5(a) and the point a = l80° in the a plot. 

For K = 2 and 3, the first term becomes dominant very early in the 
program, reflecting the constant value of V for K = 2, the linearly 
decreasing behavior of V for K = 3, and the thrust direction a = l80° 
of Vuy for both values of K. 

It was previously noted that the second and third terms of 
equation (35) 

( S ■ sin Lq — ) nu 3 X u v + 2 — u 3 • X u v 

\ Ro9s V 0 ^Rc / V 0 

tend to cancel for n = -1. These terms are due to SWp X V and 2Wj a X V 
and are directed perpendicular to uy. If n = -1 and the initial 
magnitude of SWr XV is larger than that of 2Wj a X V, which is usually 
the case for large values of S and V 0 / R o 0 s , then the combined terms 
are initially perpendicular to_ uy in the direction -u 3 X uy. It is 
clear that the magnitude of SWp X V disappears much more rapidly than 
that of 2Wq a X V, so that at some point the two terms will cancel each 
other after which the term 2Wj a X V will be larger than SW^ X V and 
the direction of the combined terms will shift by l80° to u 3 X Uy. In 
figure 5(a) this cancellation occurs near t = 0.6, and for K = 1 a 
very rapid and large shift in a is noted. The curves for K = 1 
follow closely the behavior of the combined second and third terms 
modified by small effects of the first and fourth terms of equation (35)* 
The first term is, of course, seen mainly in the impulse at t = 0. 

For K = 2 and 3, the first term is dominant and the curves follow 
closely the behavior of the first term modified by the combined effect 
of the second and third terms so that the cancellation effect is over- 
powered. Thus the early portions of all the curves reflect the rapid 
demise of the terms SWr X V (with t 3 for K = 1, t 8 for K = 2, and t 13 
for K = 3). After the cancellation, the curves for K = 1 reflect the 
rise of f* to the constant value of 2Wq a X V but for K = 2 and 3 
the term 2Wp a XV is overpowered by the first term. 

If n = +1, both Wj a and Wr are in the same direction ( see 
sketch (e)). In this case the two terms add and no cancellation effects 
occur. The case plotted in figure 5(b) has identical initial conditions 
to those for figure 5( a ) except that n = +1 for figure 5(b). As 
expected, the cancellation effects do not occur. 

The procedure of making f dimensionless with V^g also simplifies 
the discussion of propellant requirements. By specifying the parame- 
ter V o /r q 0 s , r , n, L q , S, and K the ratio of the characteristic 





27 


velocity AV to the initial relative velocity may be given: 

^ = _ JL Zn ?i£ = 0 t Pillil 

V Q V 0 Mo s f Jo V O 0 S 

Several plots of AV/v o versus Vo/Hd^s are gi ven in figures 6(a) and 6(b) 
for stated values of F 0 , n, L 0 , S, and K. In figure 6(b) an additional 
curve is given by a dotted line for a different value of r o but other- 
wise represents the same initial conditions and system parameters. As 
expected, the two curves do not differ except at low values of V o /R o 0 s . 

In summary, important influences on the thrust program from four 
parameters should be noted; namely, those of Yq/R^q, S, K, and n. 
Equation (35) consists of four terms, the fourth of which contributes 
the only position-dependent influence on the magnitude of the required 
thrust. For^large V o /R o 0 s the fourth term may be neglected; for 
small Vq/Rc^s the required acceleration for the maneuver is of the 
same order of magnitude as the gravity difference and the fourth term 
must be retained. 

The second term in (35) j due to SW-^ X V, is the largest term at 
the start of guidance for sufficiently large values of S and V o /Ro0 s . 

For usual values of S it is also the most rapidly disappearing term 
in the equation, and this behavior is usually seen In the early portion 
of the thrust programs . 

The parameter, K, primarily influences the final behavior of the 
thrust program through the first term (Vuy) in equation (35)- 

For n = -1 the second and third terms (SWr X V) tend to cancel 
each other, causing the large and rapid shift in thrust angle observed 
for cases having K = 1. For n = +1 the cancellation does not occur. 

Finally, it may be noted that some of these effects are due to the 
rotation of the reference frame, Wj a . If an inertial frame were used, 
the Coriolis terms of equation (22) would not appear, in which case 
there would be no cancellation effects for n = -1, and the limiting 
behavior prior to the impulse for K = 1 would be different. In general, 
there would be an increased uniformity among the thrust programs. 


CONCLUSIONS 


It has been recognized that the rendezvous problem is capable of 
general kinematic solution and the theory of proportional navigation 
for interception has been extended to rendezvous by modification of the 
end conditions. A simple navigation function has been selected and the 
resulting solutions of the rendezvous problem have been investigated. 
The investigation indicates the following conclusions: 
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1 . The extension of proportional navigation theory to produce 
rendezvous requires the additional parameters, Kq, which, in the simple 
function considered, reduce to the single parameter K whose primary 
effect is to control velocity so as to meet the rendezvous end conditions. 
This parameter also strongly influences the maneuver time. 

2. The thrust program requires variable thrust, both in magnitude 
and direction. The maximum value and the range of variation of the 
thrust depend on the initial conditions of relative motion and the 
constraint parameters, S and K. The behavior of the thrust program is 
strongly influenced by the value of S in the early portion, and by 
the value of K in the later portion near rendezvous . 

3 • An acceleration forcing function appropriate for use in automatic 
feedback guidance systems has been derived from the kinematics. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 1^, 196 1 
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NOTES ON THE SECOND TRAJECTORY CONSTRAINT EQUATION 


The function taken in equation ( 13 ) of the text as the second 
constraint equation has, for the purposes of the present report, the 
advantage of simplicity. However, it may he generalized to the polynomial 


w 

f(T) =^ ai ' 
i=i 


Ki(S-l) 


(Al) 


where the are taken as an ascending sequence of arbitrary constants . 

The conditions of equation (12) applied to this polynomial produce the 
following restrictions: 


S > 1 


Ki > 1 
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Without exploring fully the range of such functions, it is evident 
from (Al) and (A2) that the permissible sequences a^ at least include 
all sequences of any length made up of positive numbers whose sum is 1. 
For example, one such set is given by: 

f(T) = iy i T Ki(s_i) (A 3 ) 

i=i 

which might be described as an "average value" system. 

Equation (Al) indicates that a host of solutions of the rendezvous 
problem must exist within the confines of proportional navigation alone. 
The analysis is not carried farther here. 



30 


APPENDIX B 

GRAVITY CONTRIBUTION TO THE REQUIRED THRUST FOR 
THE SATELLITE RENDEZVOUS PROBLEM 


For the satellite rendezvous problem the effects of gravitation 
appear in equation (28) as the difference in_the acceleration due to 
gravity on the satellite and vehicle fg - f g . 


A central force field will "be assumed and the gravity forces will he 
expressed in the geocentric reference space described in the text. The 
positions of the satellite and vehicle from the center of the earth are 
R s and By, and the acceleration due to gravity for a central force field 
at the satellite is g s . 


fg s = -Ss^i 
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During terminal guidance the two vehicles are sufficiently close 
together to assume 
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and therefore (R v /R s ) 3 may be expanded binomially. The gravity 
difference then becomes: 
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Only the first-order terms will "be retained in the computations for 
the satellite rendezvous: 




*0 R , . _ _ _ v 

Ss R - "S - (“SxiUj. + x 2 u 2 + X3U3) 
-^0 


(BlO 


The successive terms in (B3) differ in magnitude by the factor 
Ro/r S j > which in the computations is restricted to values less than 0.02. 
This corresponds to values of initial range, R 0 , up to 100 miles for 
near-earth orbits. Greater ranges may he included in the computations 
by retaining the second-order terms in (A3) • It may be noted that the 
first term to appear due to oblateness is of order j(R E /R s ) 2 (RQ/R s )g s 
where R E is the earth’s equatorial radius and J, the oblateness 
constant, is about 10 -3 . This term will be less than 10 percent of the 
second-order term in (B3) for values of Ro/Rs °f 0.02 or somewhat 
greater, and therefore the gravity difference may be given up to second- 
order terms assuming a central force field. 


The accuracy with which the gravity difference need be computed 
depends on its magnitude relative to the magnitude of the remaining terms 
in the expression for the thrust program (eq. (28)), that is, the magnitude 
of the inertial relative acceleration due to external forces compared to 
the magnitude of the inertial relative acceleration required by the 
commanded maneuver. Three regimes may therefore be recognized: The 

commanded relative acceleration may be much greater than the gravity 
difference which may then be neglected; the commanded relative accelera- 
tion may be of the same order of magnitude as the gravity difference which 
must then be retained to first order; or the commanded maneuver may be so 
close to the free-fall motion that the commanded acceleration and the 
first-order gravity terms cancel each other, in which case the second- 
order gravity terms become the largest. 

For the problem discussed in the text (coplanar rendezvous with a 
satellite in circular orbit) , it was seen that the regime of operation in 
any case depends on the initial conditions; in particular, on the 
dimensionless parameter V o /Ro 0 s . If this parameter is large, then the 
case is in the first regime and the gravity difference may be neglected. 

If the parameter is small, then the case is in the second regime and the 
first-order gravity difference must be retained. In this case the required 
thrust per unit mass computed for the maneuver is necessarily small, that 
is, of the order of the gravity difference. The computations involved 
cases in both these regimes, and hence the gravity difference term was 
retained to first order. 


The third regime does not occur except under the fortuitous 
circumstance that the commanded relative acceleration is the same, both 
in magnitude and direction, as the free-fall acceleration. This event, 
if it should occur, will be of such short duration and the acceleration 
error sufficiently small that no significant motion errors can build up. 
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The second-order gravity difference term may also he important in the 
second regime from the point of view that a small thrusting error 
maintained over a sufficiently long maneuver time will huild up signifi- 
cant position errors. However, for the proportional navigation system 
this effect is obviated by fairly short maneuver times and by the feed- 
back nature of the system. The second-order gravity difference terms 
may therefore be neglected in general. 
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TABLE I.- MASS RATIOS (Mp/Mq) FOR VARIOUS SETS OF INITIAL CONDITIONS AND SYSTEM PARAMETERS *£ 

(a) T 0 = 3*/2, S = 5, ii — -1, Lq = rt/2 


IfvjjM 

0.005 

0.01 

0.02 

0.03 

0.04 



0.05 


\ K 
vA 

1 

2 

— 

3 

1 

2 

3 

1 

2 

3 

■ 

2 

3 

1 

2 

3 

1 

2 

3 

0.01 

.03 

.06 

.1 

0.973 

.929 

.861 

.776 

O.981 

.9W 

.895 

.827 

0.982 

•951 

.898 

.833 

0.955 

•927 

.863 

•779 

O.960 

.951 

.899 

.832 

0.958 

•953 

.903 

.838 

0.884 

.907 

.859 

.782 

O.87O 

.930 

.904 

.841 

0.858 

.931 

.908 

.847 

0.780 

.870 

.846 

•779 

0.738 

.884 

.893 

.845 

O.713 

.878 

.896 

.852 

O.655 

.819 

.823 

.770 

O.586 

.820 

.866 

.842 

0.550 

.807 

.865 

.848 

0.523 

.758 

•793 

.756 

0.435 

.743 

.828 

.828 

0.394 

.722 

-823 

.832 


(b) bJ b /c = 0.01, r G = 3*1/2, v 0 /c = 0.05, S = 5 


XnLo 
K \ 

- 3 * A 

-rt/2 

-rt /4 

-rt/18 

rt/18 

rt /4 

rt/2 

3 rt /4 

1 

0.901 

O.897 

0.916 

0.937 

O.922 

0.892 

O.87O 

0.870 

2 

.917 

• 917 

. 93 ^ 

•952 

.943 

.916 

.894 

.890 

3 

.919 

.920 

•937 

•953 

.945 

•920 

.899 

.894 


r 
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TABLE II.- COMPARISON OF MASS RATIOS FOR TWO -IMPULSE AND PROPORTIONAL 

NAVIGATION METHODS 


Initial 

conditions 

Two -impulse 
method 

Proportional 

navigation 

( S =5) 

to 

K 

P 

Vo 


T 

0stf,* 

Mf 

K 

0 s"kf } 

Mf 

c 

1 o 

c 


-k) 

radian 

Mo 

radian 

Mo 

0.005 

JL 

o 

0.01 

-1 

ir 

k.O 

O.985 

1.0 

0.5 

O.981 







2.0 

1.0 

.989 








3.0 

1-5 

•990 

.01 

2 

.01 

-1 

7 T 

15 

k.O 

•975 

1.0 

1.0 

.966 








2.0 

2.0 

•975 








3-0 

3.0 

.976 


* Computations for the two -impulse theory may be made 


dimensionless in the same maimer as the proportional 
navigation results. In that case it is again unnecessary 
to specify the orbit or a value of c. The free param- 
eter of the two -impulse method is no longer the maneuver 
time, tf, but the portion of the circular orbit that the 
satellite passes through during the maneuver, 0stf, which 
may also be regarded as a dimensionless maneuver, time. 




TABLE III.- EFFECTS OF POSITION ON TEE THRUST PROGRAM FOR 
LARGE V o /Rq 0 s j S = 5, K = 1 

Initial conditions: n = -1 Rc ^ s = O.Ol Ln = — 


Dimen s i onle s s 
thrust -to -mass 
ratio 


Thrust direction 
a, deg 


r o = 0 

r o = T 

r 0 = 0 

(a) 

*-1 

0 

”? 11 
ro JH: 

5.166 

5.194 

21.70 

21.58 

2.898 

2.913 

44.24 

43.75 

1.797 

1.798 

61.36 

60.56 

1.110 

1.103 

73.32 

72.16 

.6379 

.6280 

81.09 

79.51 

.3092 

.2992 

85.70 

83.16 

.0903 

.0820 

87.79 

80.75 


•3 .01+28 
.2 .1117 

.1 .1378 

. f. 1^25 


.0505 

.1165 

. 11+02 

.11+25 

(e) 


270.84 

270.34 

270.19 

'270 

180 


278.88 

272.67 

271.16 

270 

180 


= 0.7795 


= 0.7794 


Impulse at t 
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(a) K = 1 


Figure 2 
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System parameters 
S = 5 
K = I 


Mass ratio; 


n = -l 


Initial conditions 


Figure 4 


Mf 


= .9666 for r o = 0 


M f 3 7T 

= .8840 for r 0 = j2 T“ 


V 0 


R. 


= .01 


Ln " ‘ 


M 


(a) 


Vr 


Rq9 £ 


= 0.5 


Effect of position on the thrust program. 
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(a) n = -1 

Figure 5.- Effects of n and K on the thrust program, S 
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Figure Concluded 














